Non-universal tunneling resistance at the quantum critical point of mesoscopic 

SQUIDs array 
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We calculate the tunnelling resistance at the quantum critical point of a mesoscopic SQUIDs array 
in the presence of magnetic flux. We find the analytical relation between the magnetic flux induced 
dissipation strength and the Luttinger liquid parameter of the system. While the experimental 
finding for the system is around 40 — 50 mK, we find the behavior of the system even at lower 
temperatures through the analysis of renormalization group. Apart from the length scale dependent 
superconductor-insulator transition, we also predict the evidence of length scale independent metallic 
state. This study also emphasizes the importance of Co-tunelling effect. 



Introduction: Josephson junction arrays have attracted 
considerable interest in the recent years owing to their in- 
teresting physical properties. Currently such arrays can 
be fabricated in restricted geometries both in one and 
two dimensions [H4loj]. There are a few experimental 
findings in the field of mesoscopic SQUIDs array. The 
authors of Rcf. (0 ) and Ref. (@) have fabricated 
the arrays of SQUIDs junctions of different numbers in 
a single chip. They have studied the current voltage 
characteristics of the mesoscopic SQUIDs array in the 
presence of a magnetic field. They have used an exter- 
nal magnetic field (B) that tunes the effective Josephson 
coupling (Ej) between the nearest neighbor supercon- 
ductors by the following relations Ej — Ej \cos(^)\, 
where $ is the external flux and $o is the flux quantum. 
They have found magnetic field induced superconductor- 
insulator (SI) quantum phase transition at around 40—50 
mK. When the applied magnetic flux is less than the criti- 
cal value the system shows constant saturated resistance 
that may arise from the source of quantum phase slip 
centers (QPS) [4|. When the magnetic flux exceeds the 
critical value the system turns into the insulating phase 
due to the flux induced Coulomb blockade of Cooper pair 
tunneling. In this letter, we study the behavior of the 
system at very low temperature where there is no ex- 
perimental findings, i.e., around few mili Kelvin 0, [![. 
The experimental result shows length scale dependent 
superconductor-insulator (SI) transition at the quantum 
resistance but this is not the whole picture, we find the 
length scale independent metallic phase. We calculate 
the tunneling resistance at the quantum critical point 
and also show explicitly the importance of Co-tunneling 
effect to get the correct physical behavior of the system. 
Analytical relation between the flux induced dis- 
sipative strength and the Luttinger liquid param- 
eter: Here we derive analytical expression of magnetic 
flux induced dissipative strength (a) in terms of the in- 
teractions of the system. At first we derive the dissipative 
action/partition function of a quantum impurity system. 
We will see that the analytical structure of this dissi- 



pative action is identical with the mesoscopic SQUIDs 
array. 

Here we consider that the impurity is present at 
the origin where the fermions scatter from the left 
to the right and vice versa. The Hamiltonian 
describing this process is Hi = V~o(R^ (0)L(0) + 
h.c) — Vq J dx8(x)cos9(x). The total Hamiltonian 
of the system H = H + Vq J dx6(x)cos6(x,T) 
H o = ~h I uK(d x 6(x,T)) 2 + j^(d x (/)(x,T)) 2 , corre- 
sponding Lagrangian of the system is 



L = 



1 



-(d T 4>{x,T)) 2 + M(<9</>(a;,T)) 2 



2-kK J u 

+V J dx5(x)cos(d(x,T)) = L + L x (1) 

where Lq and L\ are the non-interacting and the interact- 
ing part of the Lagrangian and K is the Luttinger liquid 
parameter of the system. The only non-linear term in 
this Lagrangian is expressed by the field 9(x = 0). We 
would like to express the action of the system as an effec- 
tive action by integrating the field 8(x ^ 0). Therefore 
one may consider 6{x ^ 0) as a heat bath, which yields 
the source of dissipation in the system. The constraint 
condition for the integration is 0(t) = 9(x — 0,r). We 
can write the partition function. 

Z = I D6(x,T)e~ fo LdT 

D6{x, t)D0(t)S(0(t) - 9(0, r))e~ So LdT (2) 

Here we use the standard trick of introducing 
the Lagrange multiplier with auxiliary field A(t). 

Z = J D0(x, t) J D8(t) J DX( T )e" /<f (io+Ll)dr 

g i JgdrX(r) (0(O,t)-0(t)) 

Z = J D0(r)e~ So Lldr J D\{ T )e- lX ^ e ^ 
fD6(x,T)eSo(- L ° +iX W^» dT 

The Fourier transform of the first term of Eq.(2) is 
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-9(q,iu n )6(-q,-iu} n ) 



At first we would like to calculate the inte- 
gral: Jq dr[L — i\(r)6(0, r)], we can write this 

2 _i_ 2 2 

term as Y, q XL„ ""2^ ? i^n)9(-q, -w n ) 
-■^■(\(iu> n )9(-q,-iu n ) + X(-iuj n )0(q,ioj n ). This 
integral appears in the integral 9{x,t). This integral 
is quadratic in 9. Now we would like to perform the 
Gaussian integration by completing the square. We can 
write the result as =j- V . , „ — 2 kKv 2 , . 
In the infinite length limit one can write, 

1 \^ -kKv f dq -kKv ttK m^ttt 

VLl^q ZH? + uV - J — .. 2 , „.2„2 - — ■ «OW 
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we would like to append this result of integration in 
the second integral of Z, i.e., the integral over A. One 
can write the integrand as Yl iUn (— j£-\(iuj n )\(— iuj n ) 
+ ^(\(iu> n )9(—q, —iuj n ) + X(—ioj n )9(q, iu) n ). This integral 
is again the quadratic integral of A, therefore the Gaus- 
sian integral can be performed by completing the square. 
We get after the integration J2 iul ^9(iuj n )9(—iuj n ). 
From these analytical expression, we obtain the effect of 
bath on 9(t). The appearance of the factor uj n signifies 
the dissipation. Therefore the effective action reduces to 

S = V —9(iuj n )9(-iuj n ) + [ dxV cos9(T) (3) 

. 7TA J 



The above action implies that a single particle moving in 
the potential Vqcos9(t) subject to dissipation with fric- 
tion constant , — e?. 

Now we calculate the dissipative action of mesoscopic 
SQUIDs array. We have already proved in Ref. [4j that 
the strong coupling phase of the system is consistent with 
the experimental findings. Here we calculate the effective 
partition function of our system in the strong coupling 
phase. Our starting point is the Calderia-Legget [Tl| for- 
malism. Following reference we write the action as 



Si — So 



— Y,u m \9 m \ 2 - 



(4) 



Here, S\ is the standard action for the system with 
tiled wash-board potential 12- 3] to describe the dissi- 
pative physics for low dimensional superconducting tun- 
nel junctions, So is the action for non-dissipative part, 
a' = ^ cos |^| (the extra cosine factor which we con- 
sider in a' is entirely new in the literature to probe the 
effect of an external magnetic flux and is also consistent 
physically), the Matsubara frequency u) m = ^wfn and 
Rq (= 6.45fc£l) is the quantum resistance and R s is the 
tunnel junction resistance, (3 is the inverse temperature. 
In the strong potential, tunneling between the minima of 
the potential is very small. In the imaginary time path 
integral formalism, tunneling effect in the strong coupling 
limit can be described in terms of instanton physics. In 
this formalism, it is convenient to characterize the profile 



of 9 in terms of its time derivative, 
d9{r) 



dr 



\ eih(r 



(5) 



where h(r — Ti) is the time derivative at time r of one in- 
stanton configuration. Ti is the location of the i-th instan- 
ton, 6j = 1 and — 1 is the topological charge of instanton 
and anti-instanton respectively. Integrating the function 
over h from —00 to 00, J_ drh(r) — 8(00) — 8(— 00) = 
2tt. It is well known that the instanton (anti-instanton) 
is almost universally constant except for a very small re- 
gion of time variation. In the QPS process the amplitude 
of the superconducting order parameter is zero only in a 
very small region of space as a function of time and the 
phase changes by ±2tt. So our system reduce to a neu- 
tral system consisting of equal number of instanton and 
anti-instanton. One can find the expression for 9(oj), af- 
ter the Fourier transform to the both sides of Eq. 5 
which yields 9(ui) = — ^ i eih(iu>)e luJTl . Now we substi- 
tute this expression for 9(ui) in the second term of Eq.4 
and finally we get this term as yv F(ji — T^e^ej, where 

Fin-Ti) = ?f E m ra eMTl ~ Tj) - Wn-Tj). We ob- 

tain this expression for very small values of lu ( — > 0). So 
F(ri — Tj) effectively represents the Coulomb interaction 
between the instanton and anti-instanton. This term is 
the main source of dissipation physics of the system. Fol- 
lowing the standard prescription of imaginary time path 
integral formalism, we can write the partition function of 
the system as 0, 0, 
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(6) 



We would like to express the partion function in terms of 
integration over auxiliary field, q(r). After some exten- 
sive analytical calculations, we get 



Z = / Dq( T )e 



(~XL„ lSr<?(^)?(-^)) + (2z Xf drcosq(r)) 

(7) 

Thus by comparing the first term of the action of Eq. 
(3) and the first term of exponential of Eq. (17), we con- 
clude that the dissipative strength a and the Luttinger 
liquid parameter of the system are related by the rela- 
tion, K — 4a. 

Quantum field theoretical study of model Hamil- 
tonian of the system and explicit derivation of 
dissipative strength: 

In our previous study, we have shown explicitly that the 
mesoscopic SQUIDs array is equivalent to the array of 
superconducting quantum dots (SQD) with modulated 
Josephson coupling. We first write the model Hamilto- 
nian of SQD with nearest neighbor (NN) Josephson cou- 
pling and also with the presence of the on-site and NN 



3 



charging energy between SQD, 

H = Hji + H EC0 + H EC i- 



(8) 



Now we would like to recast our basic Hamiltonians in 
the spin language. This is valid when Eco >> Ej\. 
It is also observed from the experiments that the quan- 
tum critical point exists for larger values of the magnetic 
field, when the magnetic field induced Coulomb block- 
ade phase is more prominent than the Ej induced SC 
phase. Thus our theoretical model is consistent with the 
experimental findings. During this mapping process we 
follow Ref. ([1% Hjt = -2E J1 J2 i (SjSi +1 - + h.c), 

Hecq — EcoJ2i^i Z ■ Heci = 4EziJ2i $i Z ^i+i Z , 
Eji = E J1Q \cos(^)\ . At the Coulomb blocked regime, 
the higher order expansion leads to the virtual state with 
energies exceeding Eco- I n this second order process, the 
effective Hamiltonian reduces to the subspace of charges 
and 2, and takes the form 0-0], 



3E, 



71 



4£ , C o 



$> Z S, z 



E, 



71 



.(Si+2 Si 



h.c). 



corrections 
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Eco 

(9) 

Hci become 
) J2iSi Z S i+1 Z One can 
to as spinless fermions 
of Jordan- Wigner 



With this 

Heci - _ {±Ezi te^ 
express spin chain systems 
systems through the application 
transformation. We have transformed all Hamiltonians 
in spinless fermions which we have not present in this 
letter. In order to study the continuum field theory 
of these Hamiltonians, we recast the spinless fermions 
operators in terms of field operators by a relation fljj ]. 



ip(x) = [e lkpX tp R (x) + 



^—ikpx 



iPl{x)] where ipR^ 



and iPl(x) describe the second-quantized fields of right- 
and the left-moving fermions respectively. We would 
like to express the fermionic fields in terms of bosonic 
field by the relation ip r (x) = 4|= e~ l ~ 
r is denoting the chirality of the fermionic fields, right 
(1) or left movers (-1). The operators U r preserve the 
anti-commutivity of fermionic fields, (j) field corresponds 
to the quantum fluctuations (bosonic) of spin and 9 is 
the dual field of cj>. They are related by the relations 
4>b = 9 — <fi and 4>l — 9 + <fi. The total Hamiltonian 
is 



H 



Ho 



AE 



Z12 



{2ita) 



dx : cos(4v K <p{x)) : 



Eco 



(d x <j)(x)) dx 



(10) 



Ho is the non-interacting part of the Hamiltonian, The 
dissipative strength (ai) of this system in the absence of 
Co-tunneling effect is 



2E 



a i 



71 



2E. n 



IQEzi 



(11) 



If we consider the total effect of Co-tunneling process, 
the system reduces to Heisenberg spin chain with NN and 
NNN interactions. In this limit the dissipative strength 
(tt2) of the system is 



1 



2E 



Q'2 



71 



4 V E.n + £(4Ezi 



373.71 2 

4B C0 



(12) 



We calculate the dissipation strength by calculating K for 
both cases and then we use the relation K — 4a. This is 
the first analytical derivation of flux induced dissipation 
strength in terms of the interactions of the system. We 
consider these two processes to emphasis the importance 
of Co-tunneling effect for this system. 
Physical Analysis of Renormalization Group 
Equation and Calculation of Tunneling Resis- 
tance at The Quantum Critical Point: 
The RG equation based on the Eq. (7) is, 



dlnb 



(l-a')z 



(13) 



Following Ref._4|, we can write fugacity depends on 
length scale and temperature as, z(L) cx L l ~ a , z(T) cx 
rpa - l physical explanations based on this RG 

equation are in order. 

In our study, the resistance is evolving due to dissipation 
effect at very low temperature (few mili Kelvin, less than 
the superconducting Coulomb blocked temperature). Ac- 
cording to our calculations, for large dissipation (a' > 1), 
R(T) oc RqT^ 1 and fix > 0. Therefore at very low tem- 
perature, the system shows stable SC behaviour (region 
B of Fig. 1). and the system shows no more saturated 
resistance behaviour (region D of Fig. 1). When a' < 1, 
the resistance of system R(T) cx RqT~^ 2 and /3 2 > 0. 
So at very low temperature, the resistance of the system 
shows Kondo-like divergence behavior (region A of Fig. 
1). These stable SC and Kondo behaviour are unseen 
in the experimental findings because they have measured 
the zero bias resistance at around 40 — 50 mK. Accord- 
ing to our calculations, for large dissipation {a' > 1), 
R(T) oc RqL^ 11 and 71 > 0. Therefore the longer array 
system shows the less resistive state than shorter array 
in the superconducting phase. When a' < 1, the resis- 
tance of system R{T) oc RqL 12 , where 72 > (/3i , /3 2 , 71 
and 72 are independent numbers). So the resistance in 
the insulating state is larger for longer array system than 
shorter one. We find the dual behavior of the resistance 
for lower and higher values of magnetic field. 
From the knowledge of LL physics, we know that for a 
particular range of K , there is a metallic state in between 
the insulating phase and the superconducting phase of 
the system. The range of K depends on the nature of 
the system. Therefore we conclude that there is also a 
metallic state for a particular range of magnetic flux in- 
duced dissipation strength. This prediction is unnoticed 
in the experimental findings. At very low temperature 
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around a ~ 1 (region E of Fig. 1), the quantum phase 
slip centres poliferate and as a result a screening appears 
in the system. A detailed explanation of QPS for this 
type of system is discussed in Ref. 0] ■ Physical explana- 
tion is as follows. 

Here we consider two QPS with co-ordinates (x%, T\) and 
(^2,72), we assume that the cores of the QPS centres 
donot overlap, i.e., \x% — X\\ > xq(= £) and |t2 — ti| > 
tq(= 4-). Where £ and A are the coherence length 
and SC orderparameter of the system. We have already 
proven that the topological charges interact with each 
other logarithmically. Therefore we can write the action 
of the system as Sqps = 2S core — [ie\e2ln\x\ — x%\ ■ 
QPS with opposite topological charge attract each other 
and same charge repel with each other. Suppose we 
consider a gas of n QPS and assume that QPS cores 
do not overlap. When a current is passing through the 
system, we can write the total action of the system as 



R(T) 




°QPS 



Sint [12|- Where, 



FIG. 1: Schematic phase diagram showing the variation of re- 
sistance with temperature based on the renormalization group 
study. Region A is the Kondo like behaviour of the system, 
region B is superconducting phase, region D is the saturated 
finite resistance and region E is the metallic phase. L2 > L\ 



(14) 



Where is the distance between the two QPS at the 
site i and j and / is the current passing through the 
system. The grand partition function of the system is 
represented as a sum of all topological charges and the 
analytical form is equivalent to Eq.6, z ~ e eors . In 
the absence of current, Eq. (14) define the model for a 2- 
dimensional Coulomb gas interacting logarithmically. We 
consider the very low temperature limit (T — > 0). The 
RG equations of this 2-dimensional Coulomb gas are the 
following |# = — 47r 2 ^i 2 z 2 , = (2 — [i)z. In our case we 



can write the second equation as 



dz 



(2-/j,)z. Fi- 



nally it yields z(T) = z(A)(|r) M . When the interaction 
coefficient is less than 2, then the RG equation diverges 
at the scale T* (= z^^A). At the scale larger than £ 
the interaction between the QPS is screened and cease 
to be logarithmic, it becomes an exponentially decaying 
function (~ e- 2 ^"*^), K is the modified Bessel 
function. 

When a' = 1, i.e., $ = ($ /tt)cos -1 (-||-), the system 
has no length scale dependence SI transition at very low 
temperature. This is the critical behavior of system for 
a specific value of magnetic field. The analytical expres- 
sion for tunneling resistance at the quantum critical point 
can be calculated by comparing the expression of oti and 
ct2 with the expression of a . Rg = Rq ( 30 6 7 ^f ^ ) and 



R 



(2) 



s - -Wi^ff + M E§f x are the tunneling resis- 
tance at the quantum critical point in the absence and 
presence of Co-tunneling effect respectively. It is clear 
from our analytical derivation that the tunneling resis- 
tance at quantum critical point is not Rq but it is rather 
(2) 

non-universal. R s is proportional to the Eco as one 



would expect from the physical criteria of the system. 
Here we discuss the importance of Co-tunelling effect ex- 
plicitly. In the SC phase when a > 1, from the analysis 
of ai, we get the condition ~ 16 7 f zl > 30E J10 \cos{^)\. 
This condition is unphysical, because the all coupling 
constants are repulsive. In presence of Co-tunneling pro- 
cess, from the analysis 012, we achieve the condition of 
SC E J}0 \cos(^)\(^E J10 \co S (^)\-30) > 2* This 
condition is physically reliable . Similarly one can do the 
analysis for the insulating phase for both absence and 
presence of co-tunneling effect. 

Conclusions: we have found the non-universal tunnel- 
ing resistance at the quantum critical point of mesoscopic 
SQUIDs array. We have found stable SC and Kondo 
phase at very low temperature, which is still unseen ex- 
perimentally. We have found the length scale indepen- 
dent metallic phase in the system. The importance of 
Co-tunneling effect has been studied explicitly. 
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